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and have a unique dispersion curve. The dispersion curve controls the transmission, 
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medium. We develop a non-local hydrodynamic theory in a consistent and unified 
way. The structure of the solutions in the quantum fluid is discussed. The creation 
probabilities of all quasiparticles are derived when any one of them is incident on 
the interface. The dependencies on frequency and angular are analysed and the 
backward reflection and refraction for rotons are discussed. 
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I. INTRODUCTION 

Many physical properties of continuous media at low temperatures can be described in 
terms of quasiparticles. The dispersion relation of superfluid helium is non-monotonic and 
the quasiparticles that correspond to different monotonic regions of the dispersion curve are 
called phonons, i?^ rotons, and rotons (see Fig.l). The rotons have negative group 
velocity, i.e. their momentum is directed opposite to their group velocity. The phonons and 
rotons are observed in many experiments, such as in neutron scattering in helium [l] and 
in the direct experiments [2, 3], where beams of superfluid helium quasiparticles are created 
by a heated solid. The quasiparticles propagate in the helium, interact and reflect from 
different surfaces and quantum evaporate helium atoms from the free surface. 

These processes have been investigated both experimentally and theoretically (see for 
example [4-8]). Interestingly, R~ rotons were not detected in direct experiments until 1999, 
when they were finally created by a specially constructed source [3] and observed by quan- 
tum evaporation. All the earlier attempts to create R^ , with ordinary solid heaters and 
bolometers, were unsuccessful. 

He II quasiparticles interact with the interface and can cause their transmission, reflection, 
and conversion into each other (i.e. mode change). These are the fundamental elementary 
processes that determine heat exchange between He II and a solid, and cause phenomena 
such as the Kapitza temperature jump at a solid-liquid helium interface, (see, for example, 
Ref. ['•]). The phonons are long wavelength fluctuations of density, with the same nature 
as longitudinal phonons in a solid, and in the long wavelength limit, the problem of their 
interaction with an interface with a solid, is reduced to the problem of transmission of 
waves with some constant frequency through this interface. This problem can be found in 
textbooks on acoustics. The theory can be generalized to take into account the structure of 
real solids [9] and the structure of the real interfaces (see [10-12]). 

For short-wave excitations, rotons, there are still various models and views regarding 
their nature (see for example [13], [14] and [15]), and the problem of their interaction with 
the interface is of significant interest. The first work in this direction was [ ]. However, 
the method used there did not take into account the possibility of simultaneous creation of 
phonons and R'^ and R' rotons on the interface, and it could not distinguish between the 
R'^ and R~ rotons. 
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In the current work we consider the problem of the interaction of superfluid hehum 
quasiparticles with interfaces in a consistent and unified way. We use the approach that 
describes a quantum fluid as a continuous medium at the length scales that include the 
average interatomic distance, in terms of dispersive hydrodynamics (see [17, 18]). This was 
developed by the authors in [19]. Such an approach is feasible because in a quantum fluid 
the atoms are delocalised, as the thermal de Broglie wavelength is much larger than n~^/^ 
where n is the atom number density in the liquid. 

The idea to describe superfluid helium as a continuous medium at microscopic scales has 
been successfully used for decades. Atkins [''"] used it in the 1950s to describe the mobility of 
electrons and ions in He II. In Ref. [1 nonlocal hydrodynamics was introduced to describe 
small oscillations in superfluid helium, and in [18] it was used to describe ripplon-roton 
hybridization. The density functional approach (see for example [ ]) also utilizes this idea. 
However, the theoretical justification of this approach remained on the intuitive level until 
the work [19]. 

The nonlocality in the theory allows one to describe a medium with an arbitrary nonlinear 
dispersion relation VL{k), which appears explicitly in its equations and is the only input 
"parameter" in the theory. The possibility of using nonlocal hydrodynamics to describe a 
medium with an arbitrary dispersion relation is discussed for example in Both phonons 
and rotons are considered as purely longitudinal excitations of the same nature, only with 
different wavelengths. The quasiparticles are described as wave packets in a continuous 
medium which has nonlocal bonds and nonlinear dispersion. The problem of their reflection 
or transmission through the interface is reduced to the problem of the interaction of the 
eigenmodes of the equations of the continuous medium with the interface. 

As an intermediate step to solving the problem for all the quasiparticles of superfluid he- 
lium, we considered different special cases of VL{k). First, the monotonic dispersion relation 
+ with a > 0, was analyzed in one-dimensional [23, 24] and three-dimensional 
[25] cases. In these references, the solutions of the nonlocal equations of dispersive hydrody- 
namics in the half-space were obtained, and different effects of nonlinearity of VL{k) on the 
problem of waves transmission through the interface were distinguished. 

Then in [26] the dispersion relation such that fi^(fc^) is a cubic polynomial of k"^ was 
considered. It is the simplest expression that can approximate the distinctive dispersion 
relation of superfluid helium, in both the phonon and roton regions. The problem of the 
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simultaneous creation of phonons and rotons on the interface was solved and the interaction 
of all He II quasiparticles with the interface was described in a unified way. The failures 
of attempts to detect R~ rotons, before the experiments of Ref. [.'], were explained, and 
predictions were made for new experiments on the interaction of phonons and rotons with 
a solid. However, the use of a simple cubic approximation for the dispersion relation of 
superfiuid helium restricted the results obtained in [_ ] to qualitative or semi-quantitative. 

In the current work we present the consistent solution of the problem of the interaction of 
quasiparticles with an interface for the case when their dispersion relation is arbitrary and 
non-monotonic, so that VL^{k'^) is a polynomial of arbitrarily large degree 5*. The probability 
of the creation of each quasiparticle at the interface is derived for all cases. This work 
includes and generalizes the results of works [23]- [26], and in the special cases discussed 
previously, the expressions obtained here transform into the ones obtained earlier. So all 
results are now presented in a unified way. 

The exact expressions for the quasiparticles' creation probabilities on the interface are 
derived analytically. Those expressions are valid for all energies and angles of incidence of 
the quasiparticles that are incident on the interface. They are also valid for monotonic and 
non-monotonic dispersion relation, and for interfaces with a solid or for the free surface. All 
the qualitative results of [26] with regard to superfiuid helium, including the explanation 
why R' rotons were not detected in direct experiments until [3], are confirmed. 

The dispersion relation of surface excitations of He II, ripplons, can also be investigated 
in the same framework, but is not considered in this work, and is the subject of another 
investigation [27]. 

In the next section we consider the equations of the quantum fluid in a half-space and the 
boundary conditions that apply on the interface. Section 3 is devoted to the solution of the 
equations, and in section 4 we discuss some general consequences of the solution's structure 
and of the use of boundary conditions with regard to the problem of waves interacting 
with interfaces. Those include a generalization of Snell's law and a realization of backward 
reflection and refraction for R~ rotons. 

In section 5 we consider the process of a phonon in the solid which is incident on the 
interface, and derive the reflection and transmission coefficients. In the next section we 
calculate the partial transmission coefficients for the waves created in quantum fluid (for the 
case when there are more than one), which are the creation probabilities of the corresponding 
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quasiparticles created in the process. In section 7 we derive all the creation probabilities in 
the process of a quasiparticle of the quantum fluid incident on the interface, i.e. phonon or 
i?^ roton or i?+ roton for the dispersion of superfluid helium. We discuss the peculiarities 
of the coefficients' frequency and angular dependencies and refine the curves, obtained in 
[26] for the coefficients calculated in a rough approximation. 



II. EQUATIONS AND BOUNDARY CONDITIONS FOR HE II 

IN A HALF-SPACE 

Let us consider the half-space z>0 filled by a quantum fiuid with equilibrium density Pq 
and dispersion relation Q{k). According to the approach of [ . we introduce the variables 
of continuous medium: the velocity v, and the deviations of pressure P and density p 
from their respective equilibrium values. These quantities must obey the laws of mass 
and momentum conservation; therefore the excitations of small amplitude are described by 
ordinary linearised equations of an ideal liquid, but the relation between p and P is non- 
local, defined by some difference kernel h{r). The problem in terms of pressure P can be 
expressed as a nonlocal wave equation. When solving the problem in the half-space, the 
integration domain is limited to this half-space [26] and the problem can be brought to the 
form 

AP(r,t) = y"rfVi/i(|r-ri|)P(ri,t), x, y, tG (0, cx)), (0, oo). (1) 

zi>0 

where the dots denote derivatives with respect to time. We assume that the interface is 
sharp enough to consider that the kernel h{r) is the same in the presence of the interface as 
it is in the bulk medium. 

In the infinite medium, the integration and definition domains of Eq. (1) are infinite, its 
right-hand part is a convolution, and the kernel is related to the dispersion relation of the 
bulk excitations Q{k) through its Fourier transform (see [J!)]) 

Thus the equation (1) describes a continuous medium with dispersion relation Q{k), 
which fills the half-space z >0. We do not try to find Q{k) from microscopies, but rather 
use the experimentally measured curve as the "input parameter" of the theory and thus 
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determine the kernel h{r). This approach imphcitly contains the inner structure of the fluid. 
The dispersion Q{k) can be arbitrary. In particular, we are interested in the case, when it 
approximates to the distinctive dispersion relation of superfluid helium (see Fig.l). 

The equation (1) is supplemented by the boundary conditions on the interface z = 0. As 
the local equations of continuous medium hold on both sides of the interface, the boundary 
conditions, also local, are obtained from their integral forms in the usual way, using the 
theory of a continuous medium. 

If we consider the interface with a solid, the boundary conditions will have the form 



PI = PI 

lz=-0 \z=+0 1 



(3) 



where v^ is the z-component of the velocity of the continuous medium. 

If we consider the free surface of superfluid helium, with surface tension a, then the 
pressure at the surface should be the Laplace pressure, and for small deviations of the 
surface from equilibrium position can be written as 

where ^ is the z-coordinate of the points of the surface (see for example [IS] or [28]). 

The real interface between helium and a solid or helium and vacuum is, of course, not 
infinitely sharp; for quasiparticles, the interface can only approximately be considered fiat 
since their wave-length may be comparable with the size of the microstructure of the solid 
surface. This especially concerns rotons, as their wave-length is comparable with the inter- 
atomic distances. The possible effects of taking this into account are discussed at the end 
of section 7. 



III. SOLUTION OF THE EQUATION FOR PRESSURE IN THE HALF-SPACE 

The equation (1) was solved with the help of the Wiener and Hopf method in [23, 24] 
in the one-dimensional case, for the situation when the function fi^(A;^) is an arbitrary but 
monotonic polynomial at k> 0; then generalized in [2 )] to three dimensions. In Ref. [26] 
the case of nonmonotonic Q'^{k'^) of degree 5* = 3 was analyzed in detail. The nonmonotonic 
cubic polynomial is the simplest function that can approximate the distinctive dispersion 
relation of superfiuid helium. Therefore the use of this approximation allowed the authors 
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in [26] to solve the problem of superfluid helium phonons and rotons interaction with an 
interface, and in particular to derive the reflection and transmission coefficients for them 
in terms of roots of the corresponding cubic equation. However, due to the simplicity of 
the approximation, the results obtained were limited to semi-quantitative, especially for 
energies close to and higher than maxon the energy range below energies where the third 
degree polynomial starts to deviate from the roton part of helium spectrum significantly 
(see Fig.l). 

In the current work we use the solution of (1) for D?{k'^) a polynomial of arbitrary power. 
The physical consequences of VL{k) being essentially nonlinear were analyzed in [24, 25]. 
Those include the presence in the solution of exponentially damped waves that correspond 
to complex roots of equation VL^{k) = u"^ with regard to k; complex amplitude coefficients 
of transmission and reflection, the frequency dependencies of the coefficients, the angles of 
transmission and angles of full internal reflection. The consequences of non-monotonicity of 
fl{k) were studied in detail in [26] and these include, among others, multiple critical angles 
and backward reflection for R~ rotons. Now we have accumulated enough understanding of 
the problem to show what happens in the general case. 

The Wiener and Hopf method (see for example [29], [30]) can be used to solve the Eq. 
(1) for the function Q{k) of rather general form. Its idea is to transform to new functions, 
to turn the right-hand part of (1) into a convolution and after the Fourier transform by r 
and t, to reduce the problem to the homogeneous Riemann boundary value problem (see 
for example [31]) in the plane of complex variable k^ (see appendix for more detail). The 
problem is completely determined by its "density" 

(5) 

where u is frequency, k = + Szkz is the wave vector and k,- its projection on the plane 

The solution of the Riemann problem is standard when G is differentiable and does not 
become zero on the real axis of variable kz [29]. However it can be generalised for our case 
when it has zeros in the real roots of equation 

n\e = kl + kl) = oo^ (6) 



with regard to k^, and is always zero when /Ct- = (see appendix). 
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The simplest way to bypass the first hmitation is to shift the real roots into the complex 
plane, while preserving the index of density G. The index of G is the key parameter of the 
Riemann boundary problem, which can be calculated as the difference between the number 
of plain zeros of Eq. (6) in the upper half-plane C+ of variable and the lower half-plane 
C-. In our case G{kz) is an even function and the index is zero. Different ways of shifting 
the roots lead to different boundary problems (more details below). 

The general solution is rather complicated. In order to make the inverse Fourier trans- 
forms and ultimately derive the transmission and reflection coefficients in analytic form, we 
assume that the function Q'^{k'^) is a polynomial of power S with regard to k'^, such that 
the only real zero of Q'^{k'^) is k'^ = 0, where ~ k'^. This assumption actually does not 
restrict the generality of our consideration, because any given dispersion curve, measured in 
experiment, can be approximated, on a finite interval of k, by a polynomial, with any given 
precision. The condition that the only real zero of Q'^{k'^), when k is zero, implies that there 
are no other singularities of G on the real line. 
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FIG. 1. Dispersion relation of the quasiparticles of superfluid helium Q{k). Dots are 
experimental data [ ], and the solid line shows Q{k) obtained by fitting il^(A;^)/A;^ 
as function of k"^ by a polynomial of degree 20, so Q'^{k'^) is polynomial of degree 
5" = 21. The dashed line shows the dispersion curve analysed in [20], with S = 3. 



For arbitrary degree 5", it can be shown that the Fourier transform of the solution of 
equation (1) in terms of r and t is 

U',[h-k,MK)] 



riu [kz-ki;,{uj,kr)y 



(7) 
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The product is taken over all the roots kz = kiziuj, kr) of Eq. (6) in the upper half-plane 
C+ of the complex variable kz- The real roots are assumed to be shifted from the real line, 
with the condition that there are as many roots shifted up as there are shifted down (this 
comes from the demand that index of the density is zero). The choice which of the roots 
are shifted up (or, equivalently, down) is denoted by the index /x. It determines the specific 
Riemann boundary problem, of which (7) is the solution, and determines to within its 
amplitude C^{uj., k^), because the solution of a Riemann boundary problem with zero index 
is unique to within a multiplicative constant [2')]. The linear combination of solutions 
with all possible selections yU, gives the general solution of (1) with given u and k,-. 

As Eq. (6) is a polynomial equation of degree S with real coefficients, either with regard 
to fc^ or to kl-, the full number of roots kiz in C+, after shifting, is 5*. At a; = 0, one of them 
is kz = i\kr\. The prime superscript on the product in the nominator designates that the 
factor with this root {kz — kz) is omitted from it. This is the consequence of the singularity 
of density G in ^^ = when ^,- = 0. 

In the limit of small frequencies, when the dispersion relation is almost linear VL{k) ^ ks, 
each factor in the lower product in (7), tends to the corresponding factor in the upper 
product, with the exception of the factor with kz, which is absent in the nominator, so (7) 
tends to 

P,{kz) = -\. (8) 



Suppose that kr>uj/s. Then the two roots of Eq. (6) ±a/w^/s2 _ ^2 g^j^g imaginary, and 
there is no choice of roots shifting: due to the condition kz gC+ we have kz = i\\/kf^—oj^~/s^\. 
The solution has the form (7) and is unique. It describes the low-frequency surface excitation 



of He II, ripplon Prippi{z) ~ e 



-\k,\z 



Now suppose that kr<oj/s. Then of the two roots ±a/c^^/s^ — fc^, one should be shifted 
up. The choice kz = +\kz\ gives us the solution Pout{z) ~ exp {i\kz{oj, kT-)\z), which is 
the wave traveling away from the interface, while the other choice gives the wave traveling 
towards the interface Pin{z) ~ ex^) {—i\kz{(jJ,kT-)\z). So in this case there are two linear- 
independent solutions, and the subscript takes two values /i = out, in. The two solutions 
describe low-frequency phonons traveling in the positive and negative directions of the axis 
z. 

In the general case, the solution in coordinate space, is obtained through the inverse 
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Fourier transform, by integrating over the upper half-plane C+ of kz- If we are interested in 
solutions with a given u and k^, we take the amplitude of the form C^{u',k'^) = C^5{uj — 
uj')5{k.r — k'^) and then the solution is the sum of S monochromatic waves 

where is the projection of the radius- vector to the plane of the interface {x,y), and the 
sum is taken over all the S roots kiz in C+ (after shifting). The amplitudes are the 
residues of the right-hand part of (7) in k^z- 

The velocity is obtained from (9) through the usual relation v = —VP/pg: 

^r^{r,t■,K,u) = ^ J2 k,a^,e^(^'^(-'^-)^+'^^''^--*). (10) 
Here kj = k,- -|- Gzk-iz- 

Real solutions, which correspond to quasiparticles propagating in the medium, are ob- 
tained by making wave packets of P^, i.e. taking the inverse Fourier transform of Eq. (7) 
with the amplitude C^(ci;',k'^), that is a function that is essentially different from zero only 
when the arguments are in the neighborhood of given uj and k,-. Such solutions consist of S 
wave packets, some of them traveling and some of them damped in z>0. 

If, for some given u and fc,-, there are positive roots kf^ of Eq. (6), it means that 
there are N types of traveling waves in the quantum fluid (e.g. or R~ rotons). On the 
other hand, then of 2N real roots kiz will be shifted up from the real line. Each of the 
roots kiz in the upper half-plane of kz gives a summand in (9) and (10), so the solution 
in the quantum fluid will contain traveling waves and S — N damped waves. In terms of 
the problem of quasiparticles' interaction with the interface, as will be shown below in more 
detail, this leads to the following: when any quasiparticle is incident, all the quasiparticles 
that can be created on the interface, which conserve energy and tangential component of 
momentum, will be created with corresponding non-zero probabilities. 

IV. TRANSMISSION OF WAVES THROUGH THE INTERFACE: 

GENERAL PROPERTIES 

In the next sections we consider the problem of quasiparticles' interaction with the in- 
terface, - their reflection, transmission and mode change. We treat quasiparticles as wave 
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packets, so, as the problem is linear, it is reduced to the problem of plane waves interaction 
with the interface. However, when calculating energy transmission coefficients, we have to 
take into account that energy (as any other quantity quadratic on amplitude) is carried in 
a wave packet with its group velocity. 

The solutions both in z>0 and z<0 consist of waves with given u and k,-. Therefore in 
order for the boundary conditions to be satisfied for all t in all the points of the interface 
(x, y), the frequencies u and tangential components of wave vectors k^ on both sides of the 
interface must be equal. In terms of quasiparticles this means that energy and the tangential 
component of momentum are conserved in the processes of quasiparticle destruction and 
creation at the interface. This leads to several general consequences. 

• The wave vectors of all the incident, transmitted and refiected waves lie in one plane, 
which we may denote by {y, z), and thus the problem is reduced to two-dimensional. 
The same consideration applies to the case of the free surface. 

• R~ rotons have negative group velocity dfl/dk <0, or "negative dispersion", i.e. their 
group and phase velocities are antiparallel. Therefore for them the effects of backward 
reflection and refraction are realized, which were predicted and described by Mandel- 
stam [ ] for a hypothetical (at that time) fluid with negative group velocity. If the 
direction of propagation along the plane (x, y) of all the quasiparticles with positive 
group velocities (i.e. all except R~ rotons), incident or created on the interface, is 
along the vector k,-, the rotons will propagate in the opposite direction. The effect 
is analogous to Andreev reflection of the quasiparticles on the interface of normal and 
superconductive phases of a supercunductor, with change of signs of electric charge 
and effective mass [3 ]. Two examples are shown on Fig. 2. 

• A form of Snell's law holds for all the traveling waves involved in the process - incident, 
reflected and transmitted, so that the quantity 



is the same for all i. Here index i enumerates the types of traveling waves with different 
phase velocities in the two adjacent media (e.g. helium R~^ rotons or solid's phonons); 
Si{ijj) = u/ki{ijj) is the phase velocity of wave i and 6j>0 is its angle of propagation, 
measured from the normal. Thus all the angles of incidence, reflection and transmission 



sin &i 



(11) 
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FIG. 2. Two examples of backward refraction (on the left) and reflection (on the 
right) of Br rotons. The scheme on the left shows a phonon in the solid incident 
on the interface, and on the right a phonon of the quantum fluid is incident. The 
directions of propagation of the quasiparticles (along their group velocities) are 
shown by big arrows, and the small arrows beside them are their wave vectors. The 
directions are determined by the condition that the projections of all wave vectors 
onto the plane of the interface are equal to k^. 

are related through (11), and it also determines the number of traveling waves in each 
process. 

Suppose wave i is incident on the interface. Then if for some j ^ i Eq. (11) gives 
sinGj>l, it means that the corresponding wave is damped and the quasiparticle j is 
not created. 

Let us assign the subscripts i = 1, 2, 3 to the three types of traveling waves in the 
superfluid - the phonons, rotons and i?+ rotons. So /ci,2,3 are the positive roots of Eq. 
(6) with regard to k, with the subscripts in ascending order of their absolute values 

< ki{u) < k2iu) < k-sico). (12) 

In this work we focus on the peculiarities that the special dispersion relation of superfluid 
helium stipulates with regard to the problem of quasiparticles' interaction with the interface. 
So, for simplicity, we describe the solid by a scalar model, i.e. as an isotropic continuous 
medium with equilibrium density ps and one sound velocity Sq, and thus assign the index 
i = in Eq. (11) to the solid's phonons. 

As transverse waves can be treated in the same theoretical framework, the theory of 
continuous medium, they do not present any difficulty. However, the calculations become 
much more cumbersome, while, on the whole, the situation does not change. Due to the 
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very small transmission coefficient of the solid- helium interface (see section 7), the reflection 
coefficients hardly change at all. For the transmitted waves additional critical angles appear, 
corresponding to the sound velocity of the transverse waves. Also it should be noted that 
taking into account both the longitudinal and transverse waves in the solid, allows one to 
consider the contribution of Rayleigh waves, which give contributions to the transmission 
coefficients of He II quasiparticles into the solid at fixed incidence angles. For phonons with 
linear dispersion this problem was solved in [34]. 



V. PHONON IN THE SOLID INCIDENT ON THE INTERFACE 

Let us consider the problem of transmission of a phonon in the solid, with frequency 
u, incident from the solid side, on the interface at angle Go to the normal, into superfiuid 
helium. 

In the scale of frequencies of superfiuid helium's dispersion relation, the dispersion of the 
solid can be safely regarded as linear and its sound velocity as constant Sq = const. In the 
scalar model the solution in the half-space 2; < is the superposition of the incident and 
reflected waves with wave vectors ko = u/so. The amplitude reflection coefficient, defined as 
the ratio of pressures in the reflected and incident waves, can be expressed as 

where 

cos 60 P 



(14) 

PsSo V, 

is a generalization of the wave impedence of the interface. Due to the boundary conditions 
(3), the ratio P/vz at ^ = — is equal to its value at z = +0 and therefore is fully determined 
by the solution in the quantum fluid, in 2; > 0. This is the quantity, through which all the 
properties of the quantum fluid affect the reflection coefficient r and the solution in the solid. 

Let us now construct the "out-solution" Pout in z>0, which is realized when a wave is 
incident on the interface from the solid. This implies that a wave packet of solutions Pout 
should contain only traveling wave packets that travel away from the interface, if there are 
any. This construction of the out-solution implies the correct choice fi of roots kiz in C+. 

is the number of positive solutions of Eq. (6) with regard to k"^ for the given u and kr- 
Then among the S summands of (9) there are N with real kiz and S — N with Imkiz >0. 
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Indeed, of the 2N real roots of Eq. (6) with regard to half would be shifted up from 
the real line, and thus be included in (9). As VL{k) approximates the dispersion relation of 
superfluid helium, the possible cases to be considered are = 0, 1, 2, 3. These cases might, 
in fact, include all the problems of interest for arbitrary dispersion. 

The three roots ^1,2,3 defined in (12) are all real when w G {^rot, ^max)^ where A^ot and 
Amax are the roton and maxon frequencies respectively (see Fig. 1) . At w < A^ot and u > A^ax 
we define them by continuity. Thus ki corresponds to phonons, k2 to rotons and ^3 to 
rotons in the whole frequency range. ^2,3 at uxArot and ki^2 at uj>Amax are complex, 
among the other 5 — 1 complex roots of Eq. (6) with regard to k for u in those intervals. 

Each root ki for z = 1, 2, 3 gives a pair i^/kf — k^ of roots of Eq. (6) with regard to kz- 
Of each pair either both roots are real or the two are complex conjugate. Let us denote by 
kiz for z = 1,2,3 the three roots out of the six that enter the solution Pout, and show that 
there is only one way to choose the triplet. 

Indeed, of a complex conjugate pair, only the root in C+ can enter (7), as P{z) should be 
bounded at z>0: ki^ = a/ kf — k'^ G C+. Of the real pair, one of the roots corresponds to 
a wave traveling towards the interface, and the other, to the wave traveling away from the 
interface. 

R~ rotons have negative group velocity, i.e. for them dQ/dk < 0. Therefore, a wave 
packet constructed of plane waves with wave vectors close to k2 and their z-components 
positive, will travel towards the interface, and vice versa. As the out-solution can only 
contain wave packets traveling away from the interface, we finally obtain the following rules 
for choosing kiz for i = 1, 2, 3: 

kf, = kKLo)-e, for z = 1,2,3; 

signfc,, = (-1)^+1 ifA;2>fc2; (15) 
ki z ^ C--^ if k'^ <^ k"^ . 

All the other roots of Eq. (6) with regard to kz are complex and break up into complex- 
conjugate pairs. Of each pair there is only one root that lies in C+, and those are the roots 
kiz for i = 4:, . . . , S . Then Pout solution is built in a definite way and is unique, with kz = kiz- 
It has the form (7) with the real roots shifted up defined by Eq. (15). The rule includes all 
the cases = 0, 1, 2, 3. 

The sound velocities of most metals (e.g. Cu or Au) are much greater than the phase 
velocity of helium s{k) = fl{k)/k, so at the interface with superfluid helium it usually holds 
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that kr<ko = uj/sq ^ uj/s{k) < ^1,2,3- Therefore the case = 3 is reahzed and we have 

0<k,,<{-k2,)<k^, (16) 

So the out-solution contains three travehng waves, which means that a phonon in the 
sohd, which is incident on the interface, creates in the hehum, with non-zero probabihty, a 
phonon, a R~ roton and a roton. In the general case N depends on the function Q{k), 
values of u, kr and sq. 

In order to obtain the impedance Z, we have to calculate the values of P and v^ at 2; = 0. 
From (9), P\^^q is the sum of all residues aout,i of the right-hand part of (7), as the latter all 
lie in C+, and therefore is equal to minus the residue of P{kz) at infinity. From (10), ^z\z=o 
is the sum of residues of k^P^k^)/ (pgUj) in fcj^, and therefore we get 

P\z—:q = — res Pout{kz)] 

^ fcz— !>00 (17) 

v.L=+o = ■ /es [kzPoutikz)] . 

The residues at infinity are obtained directly by expanding (7) and on substituting them 
into (17) and (14), we obtain 

ry _ Pq f^Oz 



Ps ^ 

hz + E [f^iz - kiz{uj = 0)] 



(18) 



i=2 

where k^z = ko cos Qsoi >0 is the ^-component of the wave vector of the incident wave. 

The function of ki^ in the denominator of Eq. (18) determines the influence of the non- 
linearity of the dispersion relation Q{k) on the reflection and transmission coefficients. In 
the case of linear dispersion, it is reduced to ki z = u / s -cosG,, where 6g is the transmission 
angle defined by (11). Then Z turns into the ordinary impedance Pqs/psSo, multiplied by 
the function of angles cos 60/ cos6q. 

If we introduce the notation 

ks = irkoz] 

s (19) 

kq = kiz + Y, [kiz - kiz{uj = 0)] , 

i=2 

then the expression for the amplitude refiection coefficient takes the compact form 
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The transmission coefficient is D = 1 — IrP and therefore 



D{u, kr) 



AkgKekq 



4:pqpskozRekq 

Pqkoz + Pskq 



(21) 



2 



kg ~\~ kq 



The most interesting case for us is the interface between superfluid hehum and a sohd. In 
this shown above, = 3. Due to smallness of the parameters Pq/ ps, Si/sQ <^ 1 for 

2 = 1,2,3, the transmission coefficient is obtained in the effective hmit kg^O, so D<^1. For 
5 = 2 and = 1 the expression (21) turns into the one obtained in [ ]. For S = N = 3 it 
turns into the resuh of [26]. 

Eq. (21) is vahd for any polynomial Q'^{k'^), monotonic or not, any incidence angle or 
any density ratio Ps/ Pq- The structure of D depends, through kq, on the number of traveling 
waves A^, which is in turn determined by (11). For example, if, for some set of parameters, 
uj and incidence angle, Snell's law (11) ensures that there are no traveling waves in the 
quantum fluid A^ = 0, then kq is imaginary and D = Q. If A^> 1, then the transmitted energy 
flow is divided between several traveling waves and we need to find the partial transmission 
coefficients which give the shares of the full energy of the incident wave, that is transferred 
to each of the newly created waves. 

Fig. 3 shows D{(jj, 0) for typical parameters of a hehum-solid interface Pq/ ps, s/so = 0.1 for 
the full range of frequencies < w < 17K. We see, that at frequencies at which changes, the 
structure of D changes and the curve has kinks. These and all the subsequent plots for the 
transmission and reflection coefficients are built for the polynomial VL^{k'^) of degree 5 = 21, 
obtained by ffiting the experimental data [ ] for D?{k'^)/k'^ as function of k"^ by a polynomial 
of degree S — 1. The use of different approximation polynomials, which were obtained by 
varying the sampling of the experimental data, gives the curves that are indistinguishable 
on the given scale. 

Fig. 4 shows the part of the previous graph in the range of energies (A^ot, Amax)- Com- 
paring with the results of [26] , we see that the qualitative behavior of the probability curves 
is the same, but the numerical values are smaller by about ~ 20%. 

VI. ENERGY FLOWS AND THE PARTIAL TRANSMISSION COEFFICIENTS 



As noted above, at the interface between a solid and superfluid helium the case A^ = 3 is 
usually expected to be realized. This means that the out-solution contains three traveling 
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fia/ks (K) 

FIG. 3. Transmission coefficients for a phonon in the solid, at normal incidence 
to the interface, as functions of frequency u (in Kelvin units), for typical values 
of parameters Pq/p.s = s/sq = 0.1. The full transmission coefficient D is shown by 
the thick line, and the partial coefficients Di,2,3 are marked by numbers 1, 2 and 
3 respectively. These curves, as well as all the subsequent ones, are calculated for 
polynomial Q'^{k'^) of degree 5 = 21, obtained by fitting the experimental data for 
Q^(fc^)/A;^ as function of fc^ by a polynomial of degree S—l. 

waves, with kiz, k2z and k^z (16), which correspond to the phonon, roton and roton 
respectively. In order to find the partial transmission coefficients, corresponding to those 
waves, we need to calculate the energy flows in each of them. 

Also the conflrmation of energy conservation in the interaction process at the interface, is 
not only of abstract interest. Indeed, the flrst attempt to apply the nonlocal hydrodynamic 
description of a superfluid to the problem of quasiparticles interaction with the interface, was 
probably made in work [ ], in which the interaction of rotons with the surface excitations 
of He II was studied. In that work the solution of the nonlocal wave equation in the 
half-space, analogous to Eq. (1), was sought in the form of even functions of z, while 
extending the integration limits to inflnity. There appeared no contradiction, as the reflection 
coefficient of a roton from the free surface had absolute value of 1 either way. However, when 
applied by the authors to the problem of two adjacent media, such approach gave physically 
irrelevant solutions, in which energy was not conserved to the extent of the non-linearity of 
the dispersion relation of one of the media. The reason for this is that integral equations 
cannot be ordinarily solved by even (or odd) continuation. One can see from Eq. (1), that 
even if the values of P{z) at z<0 had physical sense, they would in fact be fully determined 
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by its values at z > 0, through the integral in the right-hand part. The demands for the 
equation to hold in z < 0, and for the solution to be even (or odd), lead to non-physical 
nonlocal bonds between the points of the fluid in the physical region z>0 with the points 
in the non-physical region z <0. This leads to a violation of the energy conservation law. 
This issue led the present authors to engage in solving Eq. (1), with the integrand given 
on the half-line z > and with finite integration limits, using the Wiener and Hopf method 
[19]-[26]. We solve it without assuming that the integrand holds at z < 0. 

Therefore let us calculate the energy flows in the solid and in the quantum fluid separately. 
First, the energy flow to the interface from the solid is the difference of the flows in the 
incident and the reflected waves. If we express the amplitudes through Pq = P|^^_q, and 
use (21), for the average 2;-component of the energy flow density we obtain 

Qoz = P^Re~k,. (22) 

Now we take into account that Eq. (6) is a polynomial equation with real coefficients with 
regard to k"^. Therefore its roots break up into real ones and into complex-conjugate pairs. 
The negative real roots give imaginary ki^, and the complex-conjugate pairs kf^ = [k'^zT 
give the pairs of roots kiz in C+, that are anti complex-conjugate: kiz = ~k*z- So all the 
non-real roots kiz do not give contributions to Refcg. The same consideration applies to the 
summands kiz{uj = ^), which are all complex. Therefore kq = YliLi ^iz and both (21) can be 



simplified and for Qqz we obtain 



IP P ^ 



Qoz = lr^y2hz. (23) 



1=1 



The expression is very simple and might seem trivial, but it turns out that thinking that 
the energy flow in the i-th wave in superfluid corresponds to the summand with subscript 
i in (23), and thus is proportional to kiz, is completely wrong. Indeed, we may remember 
that k2z < 0, while each wave created in superfluid can only carry energy away from the 
interface, not towards it. Let us calculate the energy flows in each wave explicitly. 

If Ciout,i is the pressure amplitude of a traveling wave i, and its group velocity is Ui, then 
the average 2;-component of its energy flow density is 

\kikiz\ 

Qiz = 77^^ \ alut.M , i = 1, 2, 3. (24) 
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The amphtudes aout,i are calculated as the residues of Pout (7) in kiz, and group velocities 
Ui can be obtained from the following representation of the dispersion law 

Q\k) = Ak'l[{k!-kUu^ = 0)), (25) 

i=2 

where A is a constant and kiz the same roots as in (19). After some transformations and 
using the boundary conditions (3) we can obtain 



aL.«< = P„^^n|^^.|^^^|^4^L . = 1.2.3. (26) 



'"' kj \kjz kiz kj z ~\~ kiz{oj — 0) 



Now we take into account the structure of the roots kiz, i-e that they are either real, or 
imaginary, or break up into pairs related as kiz = ~kjz- Therefore for the case = 3 we 
have 

Q^^ = \P^-k^. l''^l''l''^l'\ {^,J,k} = {l,2,3}+perm., (27) 

where perm, denotes all permutations. We can see that for any i = 1,2,3 the energy flow 
Qiz is positive, but not proportional to kiz- It can be now shown easily from (27), that in 
all the cases < 3 

N 

^^Qiz=Qoz- (28) 

2=1 

So energy is conserved and the solution is consistent. 

The relative energy flows in the three waves are functions of only kiz,k2z and k^z- This 
is very important. All the complex roots kiz are obtained by constructing an approximation 
polynomial Q'^{k'^) for the experimental data for the dispersion curve, then finding all of the 
roots of Eq. (6) on the complex plane, then sorting out those of them that lie in C^. Those 
are the necessary steps in order to, for example, build the graph D{u, k^-) given by Eq. (21) 
in Fig.3, as the complex roots enter the expressions explicitly. The functions kiz{uj, kr) differ 
slightly for different approximations, with the most rough values given by the case S = 3 
(see [2(:)]). 

However, now, in order to calculate the relative energy flows in the phonon, roton 
and -R^ roton waves, based on Eq. (27), all we need is the experimental data for A;i^2,3(i^), 
and the result is not affected by the approximation polynomial. 

The partial transmission coefficients for i = 1,2,3 are Di = DQiz/Qoz- In the most 
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interesting case A^ = 3 then we obtain 



D, = D- 



ki z~\~kj z ki z~\~ fcfe 



ki z~\~kj z~\~ kj^ z kiz kj z ki z kj^ z 



{i,j, k} = {1, 2, 3} + perm. 



(29) 



The asymptotes of the frequency dependencies of the transmission coefficients, at u close 
to the roton minimum A^ot or maxon maximum Amax, are determined by the factors that 
are functions of ki^2,3z only. When u Arot + 0, the small parameter is u = u — Arot, and 
{k2z + k3z) ~ — ^ 0, therefore D2,3 = 0{u)). Likewise when u — )■ Amax — 0, we obtain 



Di,2 = 0(v/A„aa-w) (see Fig.4). 

At normal incidence kr = and fci 3^ = fci 3, while for the R~ rotons we have 
Therefore, from (29) we obtain 



(30) 



In (30) both fractions in D2 are less than unity, due to the negative dispersion of the R~ 
rotons, which implies A;2 2 < 0, as opposed to the analogous expressions for £'1,3, in which one 
of the fractions is always grater than 1. This is the reason that D2 < D13. The graphs of 
Di/D for normal transmission are shown on Figs. 3 and 4, and there it is seen clearly that 
the difference in partial transmission coefficients is significant: 



D2 < /^LS- 



(31) 




11 12 

tialkQ (K) 

FIG. 4. An enlargement of the Fig. 3 in the range uj G (A^ot, Amax)^ showing greater 
detail. Note that the ratios Di/Dj depend only on /ci^2,3- 

The graphs on Figs. 3,4 display the same qualitative behaviour of Di{u, 0) as the analo- 
gous graphs obtained in [26], and the two sets look very similar to each other, which shows 
that the preliminary results of [-u] were correct. 
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Due to the strong inequality Ssoi ^ Si, and the generahzed Snell law (11), the wave vec- 
tors of the traveling waves of superfluid helium are concentrated in a narrow cone around 
the normal: sin Bj ^ sin Qgoi, and so Qi 71/2 for i = 1,2, 3. Therefore in the zero approx- 
imation of the small parameter s/ssoi, the angular dependencies of transmission coefficients 
are reduced to the multiplier cos6o ~ 



and thus Eq. (31) holds for any angles of incidence. 

In terms of quasiparticles this result means that a phonon in the solid, incident on the 
interface with superfluid helium, creates R~ rotons with much smaller probability than it 
creates phonons and R'^ rotons in the helium (those probabilities are small themselves, too, 
due to the small coefficient of transmission of the wave at the interface). Due to the law 
of detailed balance, likewise, when an roton is incident on the interface with a solid, a 
phonon in the solid is created with much less probability than when a phonon or i?+ roton 
in the helium, is incident. 

Thus, in an experiment in which quasiparticles of superfluid helium are created by a 
solid heater and registered on a solid bolometer, R~ rotons are weakly created and weakly 
detected (where "weakly" means much smaller than helium phonons or R^ rotons), and 
the effect gets squared. This result explains why R~ rotons were not detected in direct 
experiments until the work [■>], in which both creation and detection of quasiparticles was 
carried out by different methods. There, the R~ rotons were created by collisions in the 
bulk liquid. 

VII. HE II QUASIPARTICLES INCIDENT ON THE INTERFACE 

Let us consider the problem of a helium quasiparticle of type j incident on the interface. 
It can be a phonon j = 1, R' roton j = 2 or R^ roton j = 3. First, we discuss the 
consequences of Snell's law (11) with regard to this case, assuming, as before, the set of 
inequalities 



If a wave j is incident on the interface at angle Qj to the normal, the value of k^- = kj sin Qj 
is set and the angles of propagation of other waves i ^ j are determined by Eq. (11): 




(32) 



So > Si > S2 > S3 > 0. 



(33) 
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sinGj = Si/sj sin9j. For i<i we can define six critical angles 




i > i; i,j = 0,1,2,3. 



(34) 



If the incidence angle of wave j is greater than 6^[, it means that k^- > ki{u), so fcj^ is 



imaginary, the i-th wave is damped and the corresponding quasiparticle is not created. 

The solution in 2; > 0, which corresponds to this problem, should contain one traveling 
wave j, which travels towards the interface (or, rather, the wave packet constructed of such 
waves should be traveling towards the interface; this makes a difference for R~ rotons). Such 
a solution of form (7) is easily constructed based on the out-solution Pout- 
Let us take the solution of the form (7) and use the same rules (15) of roots shifting 
from the real line as used in construction of Pout, except for the pair ±kjz, which will be 
shifted in the opposite directions. Thus the wave packet j traveling away from the interface 
is replaced by the wave packet of the same type traveling towards the interface, and the 
notation differs from Pout by the sign of kjz (only in the denominator, as the roots kiz in 
the numerator of (7) do not change sign). Let us denote this new solution as P/^\ We can 
construct such solutions, after the number of real pairs of roots ±kiz of Eq. (6): 



The amplitude of Pj„ is Cj-„ , and the amplitude of pressure in each monochromatic wave 

(i) 

in it is j, in full analogy to the notations for Pout- 

Now, for the case A^ = 3 we have four solutions of form (7), Pout and P/;J^ for j = 1, 2, 3, 
which are linearly-independent due to their structure. The full number of solutions that 
can be formed by picking different triplets out of the real six roots to be shifted into C+ is 
(g) = 20 (where (^) is the binomial coefficient). In the case 5 = A^ = 3 it was shown explicitly 
[26], that all of them can be presented as linear combinations of the four built above. The 
solutions of Eq. (1) with the corresponding kernel were sought in the form (9), and it was 
shown that the space of solutions is four-dimensional. In the general case the additional 
complex roots kiz with i > 3 should not affect this circumstance, and we assume that, in 
quantum-mechanical terms, the four solutions form the basis set for the level at given u and 
kr- The degeneracy of the level is due to the non-monotonicity of the dispersion relation 
fl{k). The actual dimension is (A^+1) and depends on u and kr- 

The basis is constructed conveniently in such a way, that when a wave i is incident, the 
solution in 2;>0 is {Pout + Pin)- 




j = l,...N. 



(35) 
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The case N — does not correspond to the considered problem. It is reahzed when 
there are no real roots ki^'- K > ki{uj) at a; < A^ot or A;^ > k^luj) at a; > A^^x- Then the 
solution Pout is exponentially damped with z. It either corresponds to the surface excitation 
of superfluid helium, ripplon, or to the complete internal reflection of a phonon in the solid 
from the interface with He II, which is usually not realized, as mentioned above, due the 
small sound velocity of superfluid helium. The are no in-solutions and the level with cu and 
kr is not degenerate. 

The case = 1 is reahzed when there is only one pair of real ki^: a) kr < ki{u!) at 
uj < Arot, b) k3{uj) > kr > k2{u) at A^^^ < a; < A^ax, c) k^ < ks{u) at a; > A^ax- The 
) and c) correspond to a phonon or i?+ roton incident on the interface when there 
are no quasiparticles of other types with the same frequency. The case b) corresponds to 
an roton incident on the interface at angles greater than the second critical ©32, so the 
phonon and R~ rotons waves have imaginary kiz and these quasiparticles cannot be created. 
This situation is qualitatively equivalent to the problem of the interaction of a quasiparticle 
with monotonic dispersion with the interface. The level is doubly degenerate, the same as 
in the usual case, when there arc just the incident and the reflected waves. The basis is 
{Pout- Pin}: with j = 1 for a) and j = 3 for b) and c). 

The case N — 2 is realized when there are two pairs of real kiz- ki{uj) < /c^ < k2{uj) at 
Arot < < Amax- It corresponds to an W' or an R~ roton incident on the interface at 
angles greater than ©3^ and ©21 respectively. Thus the phonon wave has imaginary kiz and 
is damped. The basis is {Pout, Pin\ Pin} the solution is Pout + Pin '^ith j = 2 when 
the roton is incident and j = 3 when the R^ roton. 

Finally, the case N — 3 is realized when the number of real roots k^z is a maximum: 
kr < ki{uj) at Arot < <^ < Amax- This corresponds to a quasiparticle i of He II incident on 
the interface at angles smaller than ©^ for all j>i. So all the three waves i = 1, 2, 3 in the 
helium are traveling waves. The basis is {Pout, Pin Pin ■> Pin} the solution is Pout+P^i^ 
with J = 1, 2 or 3 depending on the type of the incident wave. 

So, when a wave i is incident on the interface, the solution in the quantum fluid is 
Pout + Pin ■ Thus it consists of one incident wave i, N reflected waves and S — N damped 
waves, with two free amplitudes Cout and C^^ (the last S waves are present in both Pout and 
P^^, so their amplitudes are summed together). The ratio C^^ /Cout is obtained with the 
help of the boundary conditions. The values of pressure P/^^ and the normal component of 
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velocity vf]^^ in the ■'-solution are calculated in the same way as for Pout, and we obtain 



P ^ 



1 r h -Lr^^'k^- 



=+0 Pouj Cout + a 



in 



where 

kg ^kiz- (37) 



If Oj < G^Q, then the solution in the solid, in 2; < 0, is one transmitted traveling wave with 
the normal component of wave vector k'^^ = — fcoCOsQo < 0, where Oq is determined from 
(11), so we have 

^-kL = - — ks. (38) 



= -0 Pq(^ 

If 9j > 6^Q then the wave in the sohd is damped and k'^^ = —i\>/kf—}^\. In this case we 
define kg such that Eq. (38) is valid again: i.e. ks = —j-k'^^. Due to the strong inequality 
■So3>Sj, B^g<^l and the first variant is realized only for helium quasiparticles incident in a 
narrow cone around the normal to the interface. 

Applying the boundary conditions (3) to (36) and (38), we obtain the ratio C^^ /Cout and 
thus express the amplitudes of all the waves through the amplitude of the incident one. The 
amplitude reflection coefficient Vij is the ratio of pressures in the created wave j and the 
incident wave i: 

nki z~ ki ;;,[u! = 0) -rr fcj^ + fcjj, ks-\-kq — 2kiz (•xo,\ 
hz + h-Xuj = Q) ■ i-i k~^z ■ k, + ~kq ' ^ ' 

'''' = ~k — k irzi — ' ^ ^- 

Kjz — Kiz kg + kq 

The energy reflection coefficient Rij is the ratio of the normal components of the energy 
fluxes in the created wave packet j and the incident wave packet i. It is also the creation 
probability of quasiparticle j when quasiparticle i is incident on the interface. In other words 
Rij is the probability of mode change from i to j, so it can be also called the conversion 
coefficient. For i=j we have Rii = and for i^j we use amplitude coefficients r^j and 
group velocities (26). Most of the factors have absolute value 1 in the similar way to the 
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transition from Eq. (24) to Eq. (27), and after some simplifications we obtain 



Rj, 



ks ~\~ ka '^ki 



ks ~\~ ka 



R 



\4:ki zkj 



■ij 



{kiz kj ■ 



{kiz~\~kkz){kj z~\~kkz') 



{kiz kkz){kjz kk 



ks~\~kq {kiz~\~kj ; 



ks + kf. 



(41) 
(42) 



where {i, j, k} = {1,2,3} + perm. 



We see that Rij = Rji. 

The formulae (41) and (42) are universal, in the sense that they are applicable to all 
possible cases = 1,2,3 and at any angles of incidence, provided kiz and kjz are real 
(otherwise it either does not have sense or Rij = 0). For example, let us find the reflection 
coefficient for the R'^ roton incident at angle 03 > B32. In this case k} = {3, 1, 2}, k^z 
is real but koz, kiz and k2z are imaginary. So only R^^ has sense. The first term in (41) is 
equal to unity. The real parts of the numerator and denominator of the second fraction are 
equal to —k^z and k^z correspondingly, while the rest of the summands are imaginary (see 
transition from Eq. (22) to Eq. (23)) and give equal imaginary parts of the numerator and 
denominator. So the second term is also equal to unity and we obtain the obvious result 
i?33 = l. Likewise, the expression for Rij (42) is simplified when k^z is imaginary (the case 
N = 2) - the middle term equals 1. 

3 

If Oj < G^Q, then kg is real and the transmission coefficient is = 1 — ^ Rij (due to the 
law of detailed balance it is the same Di, as functions of the conserved quantities u and kr, 
in Eq. (29)), otherwise kg is complex and Di = 0. The energy conservation law has the form 

3 

Rij + Di = l for < N, (43) 

including the cases when some of the coefficients are equal to zero. It was verified in [26] 
explicitly in all cases for 5* = 3. 

At the interface between superfiuid helium and a solid kg <^ ki for all i, so all the transmis- 
sion coefficients (21) and (29) are small D,Di<^l and the reflection coefficients Vij (39,40) 
and Rij (41,42) are obtained effectively in the limit kg—^O. 

The frequency dependencies of Rij and their asymptotes close to Arot and A^ax are 
determined by the factors in (41,42), that are functions only of /ci,2,3^- At uj = u—Arot^+0 
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we have {k2 z + ks z) = 0{u^^'^) , and therefore from (41) and (42) we obtain 

Rn = l-D,- Oitu'/'); R23 = I - Oitu^/^); (44) 

i?22,33 = 0(00); Ru,13 = 0(cDl/2). (45) 

So the dominating processes on the interface near the roton minimum are reflection of 
phonons into phonons, conversion of R~ rotons into R^ rotons and visa versa. 

Likewise we can derive the asymptotes of Rij at u ^ — for kr<ki^2- They can be 
obtained from (44) and (45) by changing the subscripts {1,2,3} to {3, 1,2} and the small 
parameter u to Amax~<^- The graphs of Rij as functions of frequency for kr = 0, are shown 
on Fig. 5. 
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FIG. 5. Reflection and conversion coefficients Rij at normal incidence as functions 
of w, for typical parameters Pq/ps = s/sq = 0.1. The curves are denoted by their 
respective pairs of subscripts ij. 

The angular dependencies of Rij{u,Qi) = Rij{uj,kr = ki sin Qi) are nontrivial even in 
the case of the same small parameters Pq/ps and s/sq, because the angles Gj now are not 
necessarily small. All the reflection coefficients are related through the energy conservation 
law (43), and therefore can be expressed through the coefficients -Di,2,3 and -R21, 22,23,31- The 
graphs of the four reflection coefficients as functions of the corresponding angles of incidence, 
at a flxed frequency hu = kB ■ lOK and Pq/ ps = s/ So = 0.1, are shown on Fig. 6 and illustrate 
the peculiarities of angular dependencies of Rij. 

Comparing with the previous results [ . ], we see that the qualitative behavior of the 
probability curves is the same, but the numerical values differ by ~ 10%, with the exception 
of the curves for Rn and R^-^. The probability Ru is on average almost twice as large as its 
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FIG. 6. Coefficients i?2j and R13 as functions incidence angles 02 and 0i respec- 
tively, huj = kB • 10 K, pq/ps = s/so = 0.1. The curves are denoted by their respective 
pairs of subscripts ij, and the dashed line shows R22 obtained in [26]. 



preliminary value, and -R33 is almost twice as small. This difference is clearly caused by the 
roughness of the approximation for the R'^ roton branch in [26]. 

From (42) we see that Rij^ki^kj ^ for i^j. Therefore for i>j, when the incidence angle 
is close to 7r/2, and the coefficient Rij tends to zero as y^n /2 — Qi; for i<j, when the 

incidence angle is close to the critical angle BfJ, kjz^O and the coefficient Rij tends to zero 
as i^BfJ — 0j. This gives the asymptotic behaviour of -R12, 13,23 at angles of incidence close 
to 7r/2 and of -^21,31,32 at incident angles close to the corresponding critical angles {Rij and 
Rji are the same functions of u and kj. but different functions of respective incident angles 
6j and 0j). We see that at grazing incidence, the only coefficients that do not tend to zero 
are Ra and so the dominating processes are specular (mirror) reflections; in the limit kg^O 
the transmission coefficients are zero and so Ru-^l. 

When the incident angle for some quasiparticle i increases and becomes greater that one 
of the critical angles G^J {j <i), the root kj z changes from being real to imaginary. Therefore 
the structure of the coefficients Rij changes even for those of them that do not tend to zero 
or unity. Thus, if we consider an R^ roton incident, when 62 = 621, -R21 tends to zero and 
i?22,23 have kinks. 

The numerator of a factor of i?23 in (42) is /C23 = kq — k2z —ksz, equal to kiz — k2z{^ = 
0) — k3z{^ = 0) + Yld=i{^iz ~ kiz{uj = Q))- As kiz changes from being real to imaginary, the 
real part of ^23 turns to zero. In the most simple case 5" = 3, which was analyzed in [26], 
the imaginary part without kiz is 21m.k2z{^ = ^) and is of the order of the exponent of the 
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one-dimensional kernel h{z; k^), which replaces h{r) in Eq. (1) after Fourier transforms by 
X and y. This exponent determines the effective correlation length in the z-direction for 
given kr- It was shown, that there is an angle G2 = 0° a little greater than Q^i, for which 
^23 tends to zero, and this corresponds to the penetration depth of the phonon wave j/cizl^^ 
being of the order of the effective one- dimensional correlation length in the quantum fluid. 
For this angle 9° we have R2i=D2 = because 92>O2i>02O5 also -R23 = 0. Therefore 
R22 = 1 and the R~ rotons are specularly reflected. The existence of such an angle enabled 
the authors of Ref. [•']■")] to suggest an experiment for detection of R" rotons by measuring 
the negative pressure that an R~ roton beam would exert on a membrane. The absolute 
value of the pressure is greatest when the R^ rotons are specularly reflected from it, without 
mode change, and this condition is satisfied for the incidence angle equal to Gg. 

In the general case, the effective one-dimensional kernel is a sum of many exponents, 
which are now also functions of u, but the combination of kiz that enters the numerator /C23 
still plays the role of some effective exponent. So R22 will turn to unity if the penetration 
depth of the phonon wave becomes equal to the corresponding effective correlation length. 
We can see in Fig.6, that this angle B2 indeed exists for the Q{k) that well approximates the 
dispersion relation of superfluid helium, as opposed to the rough approximation with S = 3 
used in [ • ]. Its value is about 40°, as opposed to the previous result ~30°, the peak is not 
as sharp, and at larger angles R22 is quite close to unity. This means that we only need to 
provide O2 > 40° for the negative pressure of the R~ roton beam to be measurable, though 
best results would be at 92~40°. 

Comparing the graphs in Fig.6 with the analogous graphs of [ ] shows, that at angles 
less than the critical angle the difference is small, but when the incidence angle is greater 
than the critical angle, the difference increases. For comparison, the dashed line shows R22 
from [ ] . 

If we consider reflection from the free surface, then the boundary conditions (3) are 
replaced by (4). For a solution ~ exp(ik^rT- — iwt), this expression turns into P = —akl^. 
The z-component of velocity of the surface in this case is v^ = —iuj^, and therefore we have 



PL=o- (46) 



Thus the ratio v^/P at 2; = is given not by (38), but by (46). Therefore the reflection 
coefficients from the free surface of superfluid helium are obtained by formally replacing kg 
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in Eqs. (39-42) by 



The quantity kf is imaginary, so the structure of the coefficients is the same as when D = 0, 
and Ej=i^ii = 1- 

The transition layer of hehum atoms at the interface can be treated in the usual approach 
as multiple layers of continuous media with corresponding equilibrium densities and equa- 
tions of state. This turns the problem into one of passing waves through this multilayered 
structure. 

Then the solution in the outmost "layer", the bulk fluid, is composed of Pout, Pin ^ Pin 
and Pj;^\ The one thing that changes is the boundary conditions on the interface between 
this "layer" and the next, which alters the relative amplitudes of the four solutions in the 
particular solution. This is due to changes in Z (14), or equivalently, ks (19) (or kf in 
(47)), so the transmission (21) and reflection (39-42) coefficients also will change through 
the replacement of Z by some effective quantity. 

However, when considering the process of a solid's phonon incident on the interface, the 
solution in the fluid can only be Pout, as it is the unique eigensolution that describes the 
fluid's reaction to outside influence (with new effective Z). So the relative transmission 
coefficients, i.e. the ratios D2/D1 and D2/D3, which are functions only of ^1,2,3? stay exactly 
the same as derived in this paper and Eq. (29) holds exactly quantitatively. Thus the result 
on the weak creation and detection of R~ rotons also holds. 

Likewise there are no changes to the derived asymptotes and qualitative behavior of 
the reflection coefficients Rij, as well as the critical angles B^J, which depend only on the 
frequency and angular dependences of ^1,2,3 2, and thus only on the dispersion relation of the 
bulk fluid. 

The influence of the microstructures of the interface can be taken into account by in- 
troducing the effective roughness of a solid surface. Estimates were made in [10, 11] (for 
phonons) and show that this roughness does not change the results significantly. 

It can be noted, that it is also possible to take into account the layer of weakly-adsorbed 
atoms on the interface, see [1 ]. Incorporation of this analysis into our model would also 
lead to a change of the effective impedance Z (14), in the way similar to the multilayered 
problem discussed above, with analogous consequences. 
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VIII. CONCLUSION 

In the current work we solve the problem of the interaction of superfluid helium phonons 
and rotons with interfaces in a consistent and unified way. We describe the quantum fluid in 
the dispersive (nonlocal) hydrodynamics approach, in which a medium's dispersion relation 
Q{k) enters the equations explicitly and serves as the only input "parameter". We present 
the consistent solution of the problem of the interaction of quasiparticles with the interface, 
for the case when their dispersion relation is arbitrary and nonmonotonic, so that is 
a polynomial of some degree S. 

We solve the equations of nonlocal hydrodynamics in the half-space and discuss some 
general consequences of the solution's structure with regard to the problem of quasiparticles 
creation on interface. These are a generalization of Snell's law (11) with phase velocities 
as functions of frequency, the realization of backward reflection and backward refraction for 
R~ rotons, existence of multiple critical angles corresponding to each pair of modes in the 
two the adjacent media. 

The creation probabilities, at the interface, of each quasiparticle of both media, phonons 
in the solid and phonons, R~ rotons and rotons in the superfluid helium, are derived 
when any of them are incident on the interface (21), (29), (41), (42). The obtained ex- 
pressions are valid for all frequencies, below and above the roton gap, any incident angle, 
and effectively arbitrary dispersion relation of the quantum fluid. The peculiarities and 
asymptotic behaviour of the probabilities as functions of frequency and angles are analysed. 

This work includes and generalizes the results of Refs. [2-3]- [26], in which different specific 
special cases of dispersion relation were considered, and in the corresponding special cases, 
the expressions obtained here turn into the ones obtained earlier, so they are now presented 
from a single point of view. All the qualitative results of [26] with regard to superfluid helium 
are confirmed, including the explanation why R~ rotons were not detected in experiments 
until [3]. For the suggestion of new experiments on R^ rotons detection, we have adjusted 
the optimal angle of incidence. The general expressions obtained for the transmission and 
reflection coefficients allowed us to refine the dependencies obtained in ['J( ] for the superfiuid 
helium - solid interface. 

We hope that the obtained results stimulate new experiments on the interaction of helium 
phonons and rotons with interfaces. 
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APPENDIX. WIENER AND HOPE METHOD 

Let us show how the integro-differential equation (1) is brought to a Riemann boundary 
value problem (also called Riemann-Hilbert problem). The idea of this transition is called 
the Wiener and Hopf method. 

We start from Eq. (1) and make Fourier transform by time: 

AP{r,t) = -tu^ J £rih{\v~ri\)P{vi,t), x,y,te{0,oo), ze{0,oo). (48) 

21 >0 

Let us introduce two new functions P^(r,t) and P~{r,t). The first one is defined as 

P- (r, t) = P(r, t)for z>0 (49) 
and P-{r,t) = for z < 0. (50) 

Likewise we demand for the second function 

P+(r,t) = for z>0, (51) 

and at z < it will be defined so that the equation 

A(p-(r) + P+(r)) = -cj^y £rih{\ri-v\)p-{vi), (52) 
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where integration now is made over the infinite space, should hold true on z < 0. At z > 
it holds automatically. 

Now the right-hand part of Eq. (52) is a convolution, so after Fourier transform by r 

P^(k) = [ d^re-'^'P^{r). (53) 



we have 

We are interested in P^ as functions of kz, while u and = exk^+Gyky act as parameters. 
Due to the demands (49) and (51), for the functions P^{z) that grow not faster than power 
law at infinity, their Fourier images (53) P~^{kz) and P~{kz) are analytical in the upper and 
lower half-plane of the complex plane of variable k^ respectively. So the upper index indicates 
the half-plane of the complex plane fcz, in which the respective function is analytical. Eq. 
(54) holds for real fc^, so in detailed notation we have 

P-(fc.;c^,k.)^^^|^^^^^ A:, G (-00,00). (55) 
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This equality gives a linear relation on the real line between the limit values of functions 
P^{kz) and P~{kz) on it, which are analytical in the upper and lower half-planes of the 
complex variable respectively. Therefore it defines a homogeneous Riemann boundary 
value problem on the real line (see for example [29]) with "density" 

G{k,,u,k.) = ^^^§^. (56) 

Its key parameter is the density's index, which can be calculated as difference between the 
number of plain zeros of Eq. (6) in the upper half-plane C+ of kz and the lower half-plane 
C_. As the function Q'^{kz) is even, the index of G (56) is equal to zero. 

The Riemann boundary problem with density G, that is different iable on the contour (in 
our case it is the real line), does not tend to zero on it and has zero index, has a unique 
solution to within a multiplicative constant, provided that there is an additional condition 
on its asymptotic behaviour at infinity [ ] . Our problem has two complications with regard 
to this standard case. First, G has zeros in the real roots of Eq. (6), and second, it is 
unbounded at zero when kr = 0. 

The detailed solution of this problem in the case of monotonic dispersion is given in 
appendix to [21], in the one-dimensional case, when kr = 0. It was shown there, that the 
singularity of G in zero changes the structure of the solution, which is still unique if we 
demand that P{z) is bounded at z>0. As the three-dimensional solution has to turn into 
the one- dimensional one when k^- = 0, it has the same structure, notwithstanding the fact 
that the singularity takes place for a single value of the parameter fc^. 

The simplest way to bypass the first complication is to shift the real roots into the complex 
plane, while preserving the index of density G. When Q{k) is monotonic, there are two real 
roots and so two ways of shifting while preserving the index. Thus we obtain two linear- 
independent solutions, that play the same role as plane waves with positive and negative kz 
in the usual case of linear dispersion. When Q{k) in non-monotonic, the situation is more 
complex. In ["'!■] the Eq. (1), with the corresponding kernel, was solved explicitly by using 
a different method, for the special case when Q'^{k'^) is a cubic polynomial and there are 
three positive roots of Eq. (6) with regard to /c^. It was shown there, that there are four 
linear-independent solutions with given u and k,-, and the general solution is their linear 
combination. It can be also shown, that these four solutions can be obtained by the Wiener 
and Hopf method by selecting appropriately which three of the six real roots of Eq. (6) 
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with regard to kz should be shifted up (or, equivalently, down) from the real line. All the 
rest of the (g) = 20 solutions obtained this way are linear combinations of the four. As the 
structure and dimension of solutions of Eq. (1) with given u and k,- can only depend on 
the number of possible traveling waves in them, we use the same scheme of roots shifting to 
derive all the linear- independent solutions of (1) in the case of arbitrary Q{k), which leads 
us to (7). 
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